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CJ . Abstract 



We study an exactly solvable model of D{D-i) non-Abelian anyons on a 
^ ■ one-dimensional lattice with a free coupling parameter in the Hamiltonian. For 

^^ , certain values of the coupling parameter level crossings occur, which divide the 

ground-state phase diagram into four regions. We obtain explicit expressions for 
the ground-state energy in each phase, for both closed and open chain boundary 
j~^ . conditions. For the closed chain case we show that chiral phases occur which 

^D ' are characterised by non-zero ground-state momentum. 
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1 Introduction 

Many-body non-Abelian anyonic states are characterised by non-trivial transforma- 
tion upon particle interchange, a property called braiding which underpins the theory 
of topological quantum computation [24] . For consistency, braiding transformations 
must be compatible with a set of fusion rules which govern the decomposition of 
product states in a manner that preserves the required symmetries. These compat- 
ibility requirements are automatically met if the symmetry algebra has the stucture 
of a quasi-triangular Hopf algebra [18]. A familiar example of a fusion rule arises 
in the case of spin-1/2 particles, whereby the state space for two spin-1/2 particles 
decomposes into triplet and singlet sectors. The well-known Heisenberg spin chain 
is a model that assigns different energies to the triplet and singlet sectors for neigh- 
bouring spins on a one-dimensional lattice. Analogously, one-dimensional models for 



non-Abelian anyons with nearest-neighbour interactions may be formulated by as- 
signing energies to the different nearest-neighbour sectors determined by the fusion 
rules. For the case of Fibonacci anyons, a detailed description of this approach can 
be found in [29]. Along similar lines there have been several recent studies of systems 
involving interacting non-Abelian anyons in order to better understand the properties 
of collective states in many-body systems [6, 12, 13, 16, 28]. 

One of the important features of the Heisenberg model is that it admits an ex- 
act solution, as determined by Bethe in 1931 [4]. In modern approaches the exact 
solution may be obtained using the techniques of the Yang-Baxter equation the Quan- 
tum Inverse Scattering Method (QISM) [21]. Through this framework an advanced 
understanding of the Heisenberg model (and related models such as its anisotropic 
generalisation the XXZ chain) continues to be developed in areas such as thermo- 
dynamics [20], correlation functions [7, 19], and dynamics [23, 25]. Our objective is 
to adapt this general program to derive a non-Abelian anyon system which is exactly 
solvable. Once this model is obtained, the specific goal is to determine the quan- 
tum phase transitions exhibited by the model and to investigate properties of the 
ground-state phases. 

The symmetry algebra we employ is obtained through the Drinfeld double con- 
struction [10] applied to the group algebra of the dihedral group D^. Application 
of the double construction yields an algebra denoted D{D3), which is necessarily 
quasi-triangular and consequently applicable as a symmetry algebra for non-Abelian 
anyons. The algebra D{D^) has a finite number of irreducible representations where 
the dimensions are either one, two or three. Associated with a three-dimensional 
representation we present a solution of the Yang-Baxter equation, apply the Quan- 
tum Inverse Scattering Method to determine the Hamiltonian, and solve this model 
through Bethe ansatz techniques. The Hamiltonian obtained in this manner has the 
feature that it can be expressed as a linear combination of two mutually commuting 
terms. This property leads to energy level crossings and as a result there are first- 
order quantum phase transitions at which the ground-state energy has discontinuous 
first derivative. Utilising the exact Bethe ansatz solution of the Hamiltonian we study 
the ground-state phases. Our analysis begins for systems with closed boundary con- 
ditions and we find that both time-reversal invariant and non time-reversal invariant 
phases exist, which are in evidence by computing the momentum of the ground state. 
We discover an unusual property at the boundaries at which the ground-state energy 
level crossings occur, whereby the degeneracy scales exponentially with system size. 
We also study the case of open boundaries which, due to the lack of translational 
invariance, means that states are not characterised by momenta. For this case we 
still find that the ground-state phase boundaries due to level crossing are present, 
again with exponential scaling of the degeneracy. 

In section 2 we briefly recall some basic properties of the algebra D^Ds). In 
section 3 we use a solution of the Yang-Baxter equation to derive an exactly solvable 
Hamiltonian for three different types of boundary conditions, which is solved through 
Bethe ansatz methods. To gain an understanding of the solutions of the Bethe ansatz 



equations which correspond to the ground state, we undertake numerical studies 
in section 4 and make several observations about the properties of these solutions. 
Guided by these we calculate the ground-state energy and momentum in section 5, 
and summarise our findings in section 6. 

2 Preliminaries 

We make use of the following two delta functions: 

^i^fl, i = J, S' = i,^' t =3 mods, 

\o, i^j, ' [0, i ^j mod 3, 

noting that the first delta function is not restricted to integer indices and at times 
will be used where the indices refer to elements of a group. 

We define Cij G M^xdiC) to be the matrix with a one in the ith row and jth column 
and zeros elsewhere for 1 < i,j < d. We extend these matrices by considering the 
indices modulo d, and they satisfy the relation 

fCj I, j = k mod d, 
. , , , , hJ,k,l e Z. 

U, 3 ^ K mod a, 



The algebra we will employ is based upon the dihedral group of order 6, D^. This is 
the group of symmetries of a triangle, generated by two elements a and r, with the 
presentation 

D3 = {a, rlo"^ = r^ = (ttctt = e}, 

where e is the identity element of the group. Equivalently we could consider S3, the 
permutation group on 3 objects, as it is isomorphic to D3. The Drinfeld double [10] 
of D^ is the vector space: 

D{D3) = C{gh*\g,heD,} 

where {h*\h G -D3} is a basis for the dual algebra of D^. Here we follow the notational 
conventions of [8, 9]. Multiplication and comultiplication are defined by 

gihlg^hl = 5IIII1] {g,g2)hl and ^{gh*) = J^ 9ik-'hy ® gk*. 

Details about the representation theory of D^D^), regarding construction of irre- 
ducible representations and decomposition of tensor products, can be found in [9, 15]. 
We mention that there are only eight irreducible representations, two of dimension 
one, four of dimension two, and two of dimension three. We will be concerned with an 



anyonic chain for a particular three-dimensional local state space, and consequently 
we only present the associated representation of D{D^) for this case: 
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3 The i?-inatrix and integrable models 



The zero-field six-vertex model i?-matrix is given by 

r\z) 
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- wz 
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W^^ — VL 







w~^ — w 


z^^ — z 
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w^^z^^ — wz j 



where we will set w = exp(2i7r/3). A descendant of the zero-field six- vertex i?-matrix 
is the following two-parameter i?-matrix 



R{zi,Z2) = N{zi,Z2) J2 
a,i,j=l 



.b=l 



2a(h+j) 



W{zi\h)W{z:2'\h-] 



(ii+'j,i+a ® (ii+a+j,i-, (^j 



where 



W{z\l) 



z-\ 



wz — w^ 



^-~s? 



and N{zi,Z2] 



-^wzi — w'^){w — w'^Z2). 



This two-parameter i?-matrix is a special case of the Fateev-Zamolodchikov i?-matrix 
[11], which itself is a special case of the chiral Potts /?- matrix [3]. One can also show 
that it has the symmetry of D{D-i) in the sense that 



i?(zi, Z2)^{gh*) = A' {gh*)R{zi, Z2), 



where 



A^{gh*)= J2 9k* ® g{k~'hy 

k£D3 

is the opposite comultiplication and the elements gh* are evaluated in the represen- 
tation (1). Setting zi = Z2 in (2) yields an i?- matrix previously discussed in [9]. 
The above i?-matrices are connected through the L-operator 



L{z) = ^ { (t(7* ^ei,2 + w^ *e2,i) ® ej,j + z [ei,i (g) ei_i,i + 62,2 ® ej+i,i]} . 



1=1 



Altogether we have that these operators satisfy the four Yang-Baxter equations: 

ri2ix)ri3{xy)r23iy) = r23iy)ri3ixy)ri2ix), 

ruix)Li3{xy)L23{y) = L23iy)Li3{xy)ri2{x), 

Li2{xi)Li3{xiyi)R23{yi,y2) = R23iyi,y2)Li3ixiyi)Li2{xi), 

Rl2{Xi, X2)Rl3{Xiyi, X2y2)R23{yi,y2) = R23{yi,y2)Rl3{xiyi, X2y2)Rl2{xi,X2). 

These equations are identities on the tensor product of three spaces Vi V2 ® V3, the 
subscripts of the i?-niatrices and L-operators indicating the pair of spaces on which 
the objects act non-trivially. Additionally we have the properties 

R{1, 1) = n, UR{zi, Z2)U = ziZ2R{z2\ z^^) and [R{zu ^2)]* = R{z*2, zl), 

where 11 is the usual permutation operator and *, now and hereafter, denotes complex 
conjugation. These properties imply the additional equation 

^12(3^2)^13(^2^2)^23(1/1,1/2) = R23iyi,y2)Ll^ix2y2)Ll2ix2). 

Using these operators we next construct integrable chains for three different types of 
boundary conditions. 

Periodic boundary conditions 

The most common application of the QISM is for the construction of closed chains 
of C sites with periodic boundary conditions [21]. For these systems we define the 
transfer matrices 

t(^'P){z) = tro [Loc{z)...Loi{z)] and t^^.p)^^^^ ^^) ^ ^^^ [Roc{zi, Z2)...Roi{z^, ^2)] , 

where the traces are taken over the auxiliary space Vq (note that dim Vq = 2 (3) 
for t*^^'P) (t'-^'^-')). As a consequence of the Yang-Baxter equations (3) these transfer 
matrices commute for different values of the spectral parameters z, Zi, Z2, i.e. 

[t(^'^)(x),t(^'^)(y)]=0, 

[t^''^\x),t^''^\y,,y2)]=0, 

[t^''^\xuX2),t^''''\yuy2)]=0. 

These relations imply that the transfer matrices can be simultaneously diagonalised 
and have eigenstates which are independent of the spectral parameter. This fact 
will be used for the construction of integrable Hamiltonians below. In addition, the 
transfer matrices satisfy the commutation relations 



(t(^-)(x))*,t(^-)(y) 
{t^''^\x)y ,t(''^\y„y2) 



0, 
0. 



Furthermore, as a consequence of [t'^^'P) (2;)] = t^'^'P\z) for real z, where f denotes the 
hermitean conjugate, the spectra oi r'^'^\z) and (t^^'^\z)) are the same. Finally, we 
note that by construction t^'^'^^z) is a polynomial of degree C in z, while the degree 
of each variable in t^^'^^zi, Z2) is C. 

Within the family of commuting operators generated by the transfer matrix t*^^'^) (zi, Z2) 
integrable Hamiltonians with local interactions can be constructed as a consequence 
of the observation that T = t*^'^'^''(l, 1) is the unitary operator describing translations 
by one site, i.e. TOiT^^ = Oj+i for operators Oi acting nontrivially on site i only. 
T^ = 1 as a consequence of the periodic boundary conditions. From the first order in 
a Taylor series of Int'^^'^^zi, Z2) around zi = Z2 = 1 we obtain a lattice Hamiltonian 
with nearest neighbour interactions: 



n 



il ai 


'l^inl*'"'"^.--^))- 


-PiC 


-a2 


C-l 






1^1, 


(i+i) + Hci 







_d_ 

dzo 



\n{t^^^P\z^,Z2))-l32C 



21=1,2:2=1 



i=i 



(4) 

where l^i = (^2 — \ ic^ ~^ ^V^) and cti, 0:2 have to be real for "H to be hermitian. The 
local Hamiltonian is 



H 



EB-i) 



a,b=l 1=1 



aiw'-^ + a2W '* 
(w~'' — w'-) 



ea+b+l,a+b <&> Ca+l^a, 



(5) 



with subscripts in (4) denoting which spaces the operator acts upon. We fix the 
overall energy scale by setting ai = cos{6) and 0:2 = sin(^) and separate both the 
global and local Hamiltonians into 

H = cos(^) n^^^ + sin(^) ?/(2) and H = cos(^) H^'^^ + sin(^) H^^^ 

with the property that 

['H^'\ n^'^ = 0, (6) 

which follows from the transfer matrix r^'^\zi, Z2) forming a commuting family in 
both variables. From (5) we see that the local Hamiltonians satisfy 

[/J«]^ = /7«, [iJ«]* = /j(2) and HiJ«H = if(2)_ 
The generator of translations is the momentum operator 

V = -ilnT = -iln[t^^'P\l,l)] (7) 

with eigenvalues being integer multiples of l-n j C. Due to the periodicity of the model 
and the non-cocommutativity of the D{D-i) comultiphcation (i.e. (A 7^ A^), we find 



6 



that the global symmetry of the model is broken due to the H^i interaction term, 
reducing the symmetry to D^ symmetry. The global D{D^) invariance can be main- 
tained by using a modified version of the QISM which incorprates generalised versions 
of the translation and momentum operators using braiding. 

Braided closed boundary conditions 

To obtain a model with braided closed boundary conditions [22] we define 

Z = L(0) and i? = i?(0,0). 
We construct the transfer matrices 

t^^^^\z) = tro [L^c{z)...L^i{z)U,...Uc] 
and 

These transfer matrices satisfy the same commutation relations as their counterparts 
in the periodic boundary case and have analogue properties. From t^^'^\zi, Z2) we 
construct the global Hamiltonian as was done for the periodic boundary case, yielding 

£-1 

i=l 

where H is the local Hamiltonian defined by Equation (5) and 

H, = GH(^c-i)cG-' with G = t(3'^)(l, 1) = n2ii?2i...n^(£-i)i?£(£-i) . 

We have the property that 

GHi(^i^i)G~ = i/(j+i)(j+2) and GHqG = Hu 

for i = !...(£ — 2). For this construction we see that G plays the role of the translation 
operator, and that Hq commutes with all the local Hamiltonians except those which 
act on the 1st or £th sites. This allows us to define a momentum operator using 
Equation (7). 
The operators 

bj = no+i)j^o+i)i 

are precisiely the local braiding operators for the anyonic degrees of freedom which 
satisfy the braid group relations [24]. We can interpret Hq as a braided boundary 
interaction term for sites 1 and C, since the action of G~^ = 6^i]^...6^^6||~^ is to braid 
the state at site 1 through to site C, this state interacts through H(^c-i)Ci and then 
G = hib2---hc-i braids the state back to site 1. 

Unlike the periodic model case, the global Hamiltonian with braided closed bound- 
ary conditions is invariant under the action of D{D^). This is an important property 
to maintain in anyonic models as the irreducible representations of the Hopf symmetry 

7 



algebra define the global sectors (or superselection rules) which characterise the total 
system's "topological charge" [18, 29]. Although the above models differ only in the 
boundary conditions, it is not clear that they are equivalent in the thermodynamic 
limit since for the braided case the boundary interaction is non-local. Our analysis 
below will be show, however, that the ground states are equivalent for £ — )> oo. 

Open boundary conditions 

We may also preserve the D{D^) invariance by formulating the model with open 
boundary conditions. In this case the construction of integrable models can be done 
within the framework of Sklyanin's reflection algebra [26]. To this end we require the 
additional operator 

3 

L{z) = ^{ [ei,i ® ei+i,i + 62,2 ® ej_i,j] + z (w*~^ei,2 + w^''e2,i) ® Ci^i] 

i=l 

(X L~\-{cuz)-') . 
Within Sklyanin's approach we consider the transfer matrices 

t(2.°)(^) =tro [Loc{z)...Loi{z)Loiiz)...Loc{z)] , 

t^^'°\zi, Z2) = tro [Rocizi, Z2)...RoiiZi, Z2)Rioizi, Z2)...Rcoizi, Z2)] , 

corresponding to specific representations of the reflection algebra based on unit K- 

matrices describing free ends. These transfer matrices satisfy the same commutation 

relations as in the periodic case and the spectra of t^^'°\z) and (t^'^'°\z)^ coincide. 

By construction t^'^'°\z) will be a polynomial of degree at most 2£, similarly the 

degree of each variable in t^^'°\zi, Z2) will be also at most 2£. 

Again a global Hamiltonian can be obtained by Taylor expansion of t*^^'°)(zi, ^2). 

In this case t*^^'°^(l, 1) is a scalar multiple of the identity operator, from the first order 

term we obtain 

c-i 

where H is the local Hamiltonian defined by Eq. (5). As with the previous cases both 
the global and local Hamiltonians can be split into two coupled Hamiltonians with a 
coupling parameter 6. Since the boundary conditions break translational invariance 
the notion of a momentum operator is lost. 

Bethe ansatz solution 

Starting from the Yang-Baxter equation (3) it is possible to construct functional 
relations for the transfer matrices t'-^''"-* and t^^'"^\ m G {p, 6, 0} [14]. In the following 
the fusion relation 



will be used to compute the spectrum of the transfer matrix t('^'™') . Here the functions 
f^™-\z) and g^™-\z) are for the periodic, open and braided models 



g(^\z) 
g('\z) 



g^°\z) 



{w^z + lf, 

{wz-lY, 

{l + w^zf = fp\z), 

(l-njzf = [-l)'^gip)(z), 

(l+w''z''){l-wz'') 



-11 



(1 - z^) 
■ {l-w^z^){l + wz^) 



[l+w^zf", 

[l-wzf^. 



To obtain the Bethe ansatz solution of the models we proceed by following the 
methods of [8] adapted to the fact that the transfer matrices t*'^'™-'(zi, Z2) are functions 
of two variables instead of the single- variable reduction case zi = Z2. Doing this leads 
to the pivotal result 



Proposition 3.1. Let {\j yz)} he the set of eigenvalues of t^'^'"^\z) and {(f) j{z)} 
a set of monic polynomials satisfying 

\f'\z^f'\wz) = f^"^\z)<pf'\w'z) + g^-\z)cpf\z). (9) 

// a vector v satisfies 

t^^^"'\z)v = \f'\z)v, [t^^'^^^z]]* V = \^^\z)v and t^'^"'\zi, Z2)v = A^"'\zi, Z2)v 

then the general form for A^"^\zi, Z2) is 



^(-2) 



(m) 



for Zi,Z2 G M where c-J is some constant (not uniquely defined by j and k). 

It should be noted that the functions (p^ (z) are not uniquely defined. However, if 



we restrict the degree of the 



k(^)i 



based upon the degree of t^^'^'^zi, Z2) then the 



bf 



{z) have been observed to be unique. 

For ease of notation, we hereafter drop the superscript (m) on the understanding 
that all subsequent relations hold for a fixed value of m G {p, b, 0} unless noted oth- 
erwise. From the above Propositon we see that the energy of the global Hamiltonian 
for the periodic and braided closed models will be of the form 



E 



ail 



i(l)0;.(l)-i(3 + 2v^)£ -a2i 0^1(1)01(1) -i (3 + ^v^)£ 



Noting that all of the global Hamiltonians are self-adjoint and that the energies of 
each will be of a similar form, we find that energies can be expressed as 

E = aiEj + Q;2-Efc 

r ^[0Ti(l)0;.(l)-i(3 + ^v^)£] form = p,6, (iq) 

where E- = < 

I ^[i07'(l)0;(l)-|(3 + ^V3)/:] for m = o. 

Also of interest to us is the corresponding momentum of an eigenstate in the case of 
the periodic and braided closed boundary models. The momentum is given by 

P = (Pj -Pt-i \n{cjk)) mod 27r where Pj = -i ln(0j(l)). (11) 

We remark that while P must be real it is not guaranteed that the individual Pj will 
be real. 

To proceed with the solution of the eigenvalue problem we have to address two 
problems: First, we need to determine the polynomials (l)j{z) occuring in the solution 
given in Proposition 3.1. Rewriting (9) as 



(j)j{wz) 4>j{wz) 

and using the fact that this is an equality of polynomials, the residues evaluated at the 
zeroes of (j)j{wz) in the right hand expression must vanish. Therefore, characterizing 
(pjiz) by the set of its zeroes zj. we obtain the Bethe equations for each model as 

Urn {z ~ Zk)\j{uj'^z) = 

giu'^Zk) (piojZk) ' 

Each set of non-degenerate roots {zk} of these equations provides us with a polynomial 
(pj{z). 

Second, we need to determine a pairing rule that states whether the polynomials 
4>j{zi) and 4>k{z2) can be combined to eigenvalues A{zi,Z2) of the transfer matrix 
t^^\zi, Z2) according to Proposition 3.1. Only pairs {i,j) which are allowed by this 
rule will determine energy and momentum of an eigenstate of the global Hamiltonian. 



4 Observations from numerical results of models 
with few sites 

In this section we begin to study these questions numerically for models with a few 
lattice sites. In these systems we can identify the root configurations to the Bethe 
equations (12) corresponding to the ground states of the system and, even more 
important, the patterns shown by these configurations which allow for an analysis of 

10 



the spectrum in the thermodynamic hmit £ — ?■ oo. Furthermore, our finite lattice 
analysis provides us with a basis to conjecture the pairing rules mentioned above. 
For each model we consider the following ansatz 



^A^} 



\\{z- iwyjk) 



k=l 



for some non-zero dj G N, or (f)j{z) = 1 for dj = 0.. According to (10) the energy 
associated with this function is given by 



d. 



^ EfeLi T^i^ - 1(3 + ^73) /: form=p,6, 



tEtiTZi^-t(3 + ^v^)/: form = o. 

Similarly, from (11) we can express the momentum for the closed boundary conditions 
in terms of the Bethe roots, 



Pj = -i ^ In (1 - iwyjk) 



fc=i 

We comment here that the Bethe ansatz solutions below are incomplete in the 
sense that they determine the eigenvalues of the transfer matrix t^^\zi^ Z2) only up to 
the factor Cj^, as they only give constraints on the variables yjk- This does not restrict 
its use for studies of the spectrum of the Hamiltonian since the energy expression Ej 
is not dependent on the Cjk- The total momentum (11), however, does depend on Cj^. 
We have found through our numerical studies that the Cjk are always real (they may 
be either positive or negative). Therefore their contribution to (11) is zero modulo tt, 
which is sufficient to determine whether a particular state is time-reversal invariant 
(i.e. invariant under complex conjugation, up to a phase) or not. 

The Bethe equations that we present also admit spurious solutions [8], i.e. there 
are solutions which do not correspond to actual eigenvalues of the transfer matrices. 
Below, the statements concerning the properties of the Bethe roots relate to solutions 
which are not spurious solutions. 

Periodic boundary conditions 

Using our ansatz for <i)j{z) the Bethe equations become 

di 



n 



wyji - yjk \ ( iwyji - 1 



w^yji-yjkj \iui'^yji 



k=l 

for < / < dj. These are the same equations found in [8]. 

Observation 4.1. We make the following observations about the Bethe roots of the 
function (f)j{z), 

11 



1. ifyjk is a Bethe root then so is y*f,, 

2. there are C or C — 1 distinct Bethe roots, one of which may he zero. 

3. any two functions (pjiz) and (f)k{z) may pair if and only if both functions have 
the same number of non-zero and zero roots. 

Of particular importance are the possible ground states, i.e. the eigenstates of H^^^ 
(or equivalently "H*^^-*) with highest and lowest energy. Our numerical results for the 
corresponding configurations of Bethe roots for the highest energy state in a system 
with £ = 2, 3, 4 are shown in Table 1. From these data we conjecture that this state 

Table 1: Bethe roots for highest energy state of "H*-^^ 



c 


Energy 


Vjk 


MVjk) 


2 


2.31 


-0.57735 + 0.81650 
-0.57735 - 0.81650 


2.1863i 
-2.1863i 


3 


1.97 


-0.38896 + 0.76875i 
-0.38896 - 0.76875i 


-0.14902 + 2.0392i 
-0.14902 -2.0392i 


4 


3.15 


-0.74151 + 1.1157i 

-0.74151 -1.1157i 

-0.41319 + 0.62170i 

-0.41319 -0.62170i 


0.29240 + 2. 1574i 

0.29240 -2.1574i 

-0.29238 + 2. 1574i 

-0.29238- 2. 1574i 



is determined by complex conjugate pairs of roots, the ln{yjk) form so-called 2-strings 
[4, 27]. The pairs are centered on the real line and seperated by 47ri/3 up to finite- 
size effects. Our data for £ = 5 and 6 support this conjecture and we have verified 
numerically that such a configuration solves the Bethe equations for lattices with up 
to £ ^ 1000 sites. 



The Bethe roots for the lowest energy states are given in Table 2. In this table 
we have included a case where two different sets of Bethe roots give the same energy. 

We observe that all of the Bethe roots yjk are real, or equivalently, the logarithm 
of all Bethe roots lie either on the real line or are shifted by in. In the case of 
£ = mod 4 it has been found that the set of Bethe roots is invariant under inver- 
sion {yjk} ^ {^/Vjk} with £/4 positive and 3£/4 negative Bethe roots. Again we 
have obtained data for £ up to 640 sites in support of this observation. 

Braided closed boundary conditions 



Using our ansatz for 



the Bethe equations become 



n 

k=l 



wyji - yjk 



w'^Vji 



yjk 



1 - iwyji 
iw^yji + 1 
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Table 2: Bethe roots for lowest energy state of 'H^-^'^ 



c 


Energy 


Vjk 


HVjk) 


2 


-1.15 


-0.86603 


-0.14384 + i7r 


2 


-1.15 




-1.1547 


— CXO 

0.14384 + i7r 


3 


-1.73 


0.50771 
-0.77786 
-1.4619 


-0.67784 

-0.25121 + m 

0.37974 + in 


4 


-2.31 


1 

-1 
-1.7321 
-0.57735 




in 

0.54933 + in 

-0.54933 + in 



For C even these coincide with the equations for the periodic case, although the 
conditions imposed on non-spurious solutions are different: 

Observation 4.2. We make the following observations about the Bethe roots of the 
function (f)j{z), 

1- if Vjk is a Bethe root then so is y*^, 

2. there are C or C — 1 distinct Bethe roots all of which are non-zero, 

3. if \{z) is an eigenvalue oft^'^'^\z) then z^\(z^^) is an eigenvalue oft^'^'^^z), 

4- any two functions <t>j{z) and (pk{z) may pair, 

5. the energy spectrum w-^' with braided closed boundary conditions is a subset of 
the energy spectrum of 1-L^^^ with periodic boundary conditions. 

Note that while the third observation implies the fifth one it does not provide any 
insights into the degeneracies of the spectrum for the global Hamiltonian "H. The 
difference in pairing rules compared to the periodic boundary model is related to the 
change in the symmetry of the global Hamiltonian i.e. the difference between D^ and 
D[Ds) symmetry. The condition that all Bethe roots be non-zero renders some of the 
configurations found in the case of periodic boundaries spurious. The solutions to the 
Bethe equations with yjk 7^ for all k for the ground states identified above, however, 
do correspond to ground states in the presence of braided closed boundary conditions. 

Open boundary conditions 

Using the ansatz for (/)j{z) the Bethe equations become 



n 

fc=i 



wyji - yjk 



w^Vji 



yjk 
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for < / < dj. 

Observation 4.3. We make the following observations about the Bethe roots of the 
function (t)j{z), 

1. if ijjk is a Bethe root then so is y%, 

2. ifyjk is a non-zero Bethe root then so is y^^ , 

3. there are 1C or 2C — 1 distinct Bethe roots, 

4- if there are 2C Bethe roots they are all non-zero, 

5. if there are 2C — 1 Bethe roots one of them is zero, 

6. any two functions (j)j{z) and (pkiz) may pair. 

Our numerical analysis of the Bethe equations for systems with a few sites shows that 
apart from the doubled number of roots the ground states are again described by 2- 
strings and real yjk, repsectively. The existence of these solutions has been verified 
for systems with several hundred sites. 

5 Ground-state energy in the thermodynamic limit 

We now calculate the ground state energy in the thermodynamic limit for the periodic, 
braided closed boundary, and open models. We begin by considering the states with 
highest and lowest energies of the Hamiltonian 1-L^^\ From our analysis of small 
systems we know that the Bethe root configurations for the highest and lowest energy 
states for the different boundary conditions are essentially the same (apart from the 
doubling of roots in the open boundary case). 

Expanding the energy in C one has the general form 

E = £eoo + Poo + -^ X const. + o(£"^), 
where £eoo and poo are the bulk and boundary energy contributions, respectively. 

Bulk Energy of H^^^ 

Here we calculate the bulk energy for the highest and lowest energy states of H^^^ 
in the thermodynamic limit £ — )> oo. By definition this term is linear in the system 
size C and independent of boundary contributions. Below we assume C to be even, 
therefore our results apply to chains with both periodic and braided closed boundary 
conditions having the same Bethe equations. 

Motivated by our numerical results for the highest energy state of the model with a 
small even number of lattice sites, we consider a distribution of Bethe roots given by 

y2k-i = 6"''=+ 3 and t/2fc = 6""= 3 
14 



for 1 < A; < £/2 with x^ G M. This configuration of Bethe roots is the aforementioned 
2-string hypothesis. In terms of the new variables Xk the Bethe equations become 



c 

2 

n 



sinh(2^ + f ; 



sinh(^ 



-X). ^7^^ 



sinh(f + iff)sinh(f + 



lliTT -\ 

12 ^ 



sinh(f -iif)sinh(f -iff 



for < / < C/2. For £ — ?> oo the roots x^ can be described by their density p{x) 
which, after a Fourier transformation, satisfies the hnear equation 



where 



p{v) 



^<"-*)=2?^ 



-B V 



l-K 



B V 



IHtt 
^2~ 



B { v;'-^ ] Piv) 



' d , /sinh(f +t) 
^ m ' 



dx \sinh(| — t) 

Solving this we find the density 

3 1 



'dx 



sinh(t;(7r + 2it)) 
sinh(7rf) 



p{x) 



2rc cosh(3x) 



We calculate the energy to be 



Ctr 



k=l 



+ 



1 - ie^k I _ ie'k+^ 



iC 



c 



p{x) 



1 — ie^ 



+ 



1 — le ^ 3 



-£(3 + ^73) 





dx £(3 + iv3) 
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This is the bulk energy of the highest energy state of "H*-^-*. 

For the lowest energy state of "H^^' we consider chains of length £ = mod 4 where 
according to our conjecture above the root configuration is given by 



Vk+i 






1 < A; < £/4, 
1 < A; < 3£/4, 



with x^jX^. G M. We define two different density functions pa^b for the two subsets 
of Bethe roots. From the Bethe equations we can construct coupled equations which 
are linear in Fourier space. Solving these equations yields 



Pa[X) 



3y2cosh(; 



■X] 



4 7rcosh(3a;) 4 7rcosh(3a;) 
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and 



Pbix) 



3y2cosh(^ 



X] 



4 7rcosh(3a;) 4 7rcosh(3a;) 
With these densities we obtain for the bulk energy of the lowest energy state of "H*^^-* 



Cer 




1 + iwe^ I 6 



Jv^ 



Boundary Energy of ^^^^ 

Next we compute the boundary energy p^ for these states. As a consequence of 
translational invariance poo vanishes in the closed chains. We can deal with open 
boundaries as above with two modifications: due to Observation 4.3 we now have to 
consider 2£ roots which for the highest energy state of "H*^^-* amounts to parametrizing 
the roots as 



y2k-\ 



3*^fc~r 3 



and 



l/2fc 



o^fc- 



2i7r 
3 



for 1 < A; < £ with x^ € IR. In addition only even distributions of the x^ are to be 
considered (see also [2, 17]). Proceeding as above we introduce a density function 
p(x) which in Fourier space satisfies 



p[v) 



-2 B 



ITX 



B 



1 



B 



V 2i7r\ 
2'~J 



B 



llin 
^2~ 
V bin 
2'"6~ 



B 






B 



ITT 



p{v) 



To order £° this equation is solved by twice the density function found for the closed 
boundary conditions above. This reproduces the bulk energy of the highest energy 
state. The order C~^ correction of the density gives the boundary energy for this 
state 



Po 



3 2v/3 

■2 + — 



For the calculation of the boundary contribution to the lowest energy state of W^' 
we proceed in the same way. As a result we obtain 



3 2\/3 
"4 + — 
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Energy, momentum and degeneracy of the ground state 

For the highest and lowest energy states we have calculated both the bulk and bound- 
ary energy analytically for "H^^^ Our numerical studies of small systems indicate that 
the excitation spectrum is gapless. The finite size corrections (of order C~^) to the 
ground state energies can be obtained through numerical analysis of the Bethe equa- 
tions for systems with up to £ = 1000 lattice sites. For the closed models we find 
that the highest energy is 



E^ 



high 



while the lowest energy is 



1 2^3 



n 



9 



£ 



12 TT ,^ -I, 



E, 



low 



1 

2^ 



2^3 3 



C 



3 vr / ^ ix 



Similarly, we present the energies for the open model. The highest energy is 



-^high 



1 2^3 

71 



9 



£ 



3 2VS 

--H 

2 3 



12 n 
~5^2AC 



o{C 



-i\ 



while the lowest energy is 



E 



low 



1 

2^ 



2^3 
~9~ 



C 



3 2v^ 

4 + — 



3 TT 
- X 

2 24£ 



o{C 



-i> 



Note that this finite size scaling behaviour is consistent with the predictions of con- 
formal field theory [1, 5]. To identify the Virasoro algebras corresponding to these 
sectors of the model we need to include the low-energy excitations in our analysis. 
This will allow us to compute the (non-universal) Fermi-velocity needed to extract 
the central charge from the £~ ^-corrections to the ground state energies and also the 
conformal dimensions appearing in the excitation spectrum. These questions will be 
addressed in future work. 

From the previously stated properties we know that the energies given above must 
also be the highest and lowest energies of 'H*^^-'. Furthermore, based on our analysis 
of the models with few sites we expect that both the highest and lowest energies 
pair. Therefore, the ground-state energy will depend upon the values of the coupling 
parameter 9 and is specifically given by. 



low 



77 
ground state 



' ^^mEio^ + sin(^)Ei^ 
cos(^)%gl^ + sin(^)EiQ^, 
cos(^)Ej^jgl^ + sin(^)Ej^igj^, 



cos 



0< ^ < 






2' 

<7r, 

2 ' 



< 



W^low + sin(^)i5high' f < ^ <2 



7r 



where we have made use of (6). 
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Due to the manner in which all the (pki^) pair to produce the transfer matrix 
spectrum when the D{D^) symmetry is not broken, it is anticipated that the number 
of (non-spurious) solutions of the Bethe ansatz equations scales exponentially with C. 
Our numerical calculations support this view. It is clear that level crossings will occur 
when 9 = rmr /2 for m & Z yielding either ai = or 0:2 = 0. These level crossings 
at which the ground state is degenerate give first order quantum phase transitions, 
as it is straightforward to show that the energy has discontinuous first derivative. 
The ground-state energy at these transition points is determined by a single solution 
of the Bethe ansatz equations, with the second solution in (10) chosen arbitrarily. 
Hence the ground-state degeneracy will depend on the number of solutions of the 
Bethe ansatz equations, leading to an unusual exponential scaling of the degeneracy. 

While the ground state is non-degenerate for 6 ^ m7r/2, it is important to briefly 
discuss the fact that D{D^) admits two distinct one-dimensional representations [9] 
which define two different global topological charge sectors. While we cannot deter- 
mine which of these sectors the ground state is a member of generally, we can establish 
that the ground state for 7r/2 < 9 < ti and the ground state for 37r/2 < 9 < 2Ti be- 
long to the same sector. This follows from the properties i) ['H*^^-']* = 1-L^'^\ meaning 
that the ground state for 7r/2 < ^ < tt is the complex conjugate of the ground state 
for 37r/2 < 9 < 27r, and ii) complex conjugation leaves all topological charge sectors 
invariant since the representation (1) is real. To show that the ground states for 
7r/2 < 9 < Ti and 37r/2 < 9 < 271 are distinct in the case of the closed models, it 
remains to finally consider the momentum. 

Using the densities of Bethe roots calculated above we can compute the momentum 
of the ground state for the closed chain models: 

R'e(i^low) = -^^ and Re(Phigh) = —C 

We have ignored the imaginary components as they have been observed to be cancelled 
out by the imaginary component contributed by Cjk- Furthurmore we observed that 
the constants Cjk are real, this determines the momentum modulo vr: 

mod vr, < ^ < f , 

p _ ; ^^ mod TT, I < ^ < vr, 

ground state = \ q mod vr, vr < ^ < f , 

-^C mod TT, ^ < 9 <2ti. 

From the above we see that a signature of the ground-state phases for 7r/2 < 9 < -k 
and 37r/2 < ^ < 27r is that the momentum is generally non-zero modulo tt, so the 
states are not time- reversal invariant. Thus these two phases are two distinct chiral 
phases that belong to the same topological charge sector. Conversely the ground 
states for < 6* < 7r/2 and -n < 9 < 37r/2 are time- reversal invariant. 



6 Summary 

The Hamiltonian we have studied admits gapless excitations in the thermodynamic 
hmit for all values of the coupling parameter 9. At the particular values 9 = rmr/2 
for m G Z ground-state energy level crossings occur, which divide the ground-state 
phase diagram into four regions. We have computed the ground-state energy for 
each of these regions for both closed and open boundary conditions. From these 
explicit expressions it is seen that the level crossing points correspond to first order 
transitions where the first derivative of the ground-state energy is discontinuous. We 
also computed the ground-state momentum in the closed chain case, to establish that 
chiral phases exist which are characterised by non-vanishing momentum. 
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